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Kinematics of two dwarf galaxies in the NGC 6946 group 
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Abstract. We present high velocity resolution (~ 1.6 km s~^) Giant Meterwave Radio Telescope (GMRT) HI 
21 cm synthesis images of the dwarf galaxies KK98 250 and KK98 251, as well as optical broad band images of 
KK98 250. We find that, despite being faint (Mb < —14), both galaxies have regular velocity fields. In both cases 
the velocity fields are consistent with rigid body rotation. We derive rotation curves for the galaxies and fit mass 
models to them. For both galaxies, we find acceptable fits using isothermal halos while NFW halos provide a poor 
fit. Finally, we compile from literature a sample of galaxies with HI synthesis observations and I band magnitudes. 
From this sample we find, in agreement with earlier studies (which used single dish HI data), that dwarf galaxies 
tend to lie below the I band TuUy-Fisher relation defined by brighter galaxies. 
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1. Introduction 

There has been renewed interest in the kinematics of dwarf 
galaxies for two reasons. (1) There appears to be a conflict 
between hierarchical galaxy formation model predictions 
for, and observational determinations of, the dark matter 
halo density profile of galaxies. In particular, these models 
predict a cuspy central halo, which many observers find in- 
consistent with observations (e.g. IWeldrake et al. 2003|l . 
While the models predict cuspy halos for all galaxies, 
dwarf galaxies are believed to be best suited for tests of 
the models because these galaxies are in general held to be 
dominated by dark matter, even in the innermost regions. 
In larger galaxies, uncertainties in the stellar mass to light 
ratio limit one's ability to determine density profile of the 
dark matter halo. (2) It appears that faint dwarf galaxies 
deviate significantly from the TuUy-Fisher relationship as 
defined by bright galaxies. While the TF relation is gen- 
erally studied using the inclination corrected 50% profile 
width (usually obtained from single dish observations), it 
is unclear whether this is an appropriate measure to use 
for the faintest dwarf galaxies, where the rotational veloc- 
ities are comparable to the velocity dispersion. For such 
galaxies, the "asymmetric drift" corrected rotational ve- 
locity is probably a better measure. However, estimating 
this quantity requires HI synthesis imaging. Further, for 
the faint dwarf irregulars, even the inclination may be dif- 
ficult to estimate from the optical image and may only 
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be obtainable from the HI distribution and kinematics. 
Unfortunately, the number of faint dwarf galaxies with 
sufficiently high quality HI synthesis images available is 
quite small. We discuss here GMRT HI observations of 
two dwarf galaxies, KK98 250 and KK98 251, in hght of 
both these issues. In discussing the TF relation, we com- 
bine our data for KK98 250 and KK98 251 with our earlier 
data for still fainter galaxies and data for brighter galaxies 
available from the literature. 

Both KK98 250 and KK98 251 were identified 
as companions to the giant spiral galaxy NGC 6946 
by Karachentseva & Karachentsev (1998) and have also 
been observed in HI at Effelsberg (Huchtmeier et al. 1997) 
and at the DRAG (Pisano & Wilcot s 2000|l . The data pre- 
sented here are of much higher sensitivity and resolution 
than the DRAG data. Distance estimates to KK98 250 
and KK98 251 vary from 5.3 Mpc (Sharina et al. 1997, 
Karachentsev et al. 2000) to 8.2 Mpc (Sharina et al. 1997). 
Because of the proximity of the two galaxies to each other 
and to NGC 6946, on the sky as well as in the velocity, 
we feel that it is likely that all these galaxies belong to 
the same group viz. the NGC 6946 group. Hence, in this 
paper we take the mean distance to the group (5.6 Mpc, 
estimated from the brightest stars in eight members of the 
group) as the distance to both the galaxies (Huchtmeier 
et al. 2000). At this distance, the absolute blue magnitude 

for KK98 250 and KK98 251 are Mb 13.72 and -14.54 

respectively. 

The rest of the paper is divided as follows. The GMRT 
observations are detailed in Sect. |21 while the results are 
presented and discussed in Sect. 13 
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2. Observations 

2.1. Optical Observations and analysis 

CCD images of KK98 250 and KK98 251 in the Bessell I 
and V filters were obtained on 29*'' and SO*"" May 2003, 
using the HFOSC (Himalayan Faint Object Spectrograph 
Camera), at the 2.0 m Himalayan Chandra Telescope. The 
camera has a field of view of lO'xlO', with a scale of 
0.3 /pixel. Both the galaxies were covered in the same 
pointing. The total exposure time on the target was 40 
min in I and 50 min in the V band. The FWHM seeing of 
the co-added images was ^1.7 . Since the first night was 
not photometric, standard fields (from 
were observed only on the second night 

Debiasing, flat-fielding and cosmic ray filtering were 
done in the usual manner, using standard IRAF routines. 
The exposures taken on the 30**^ were calibrated and then 
added after alignment. This combined frame was used to 
calibrate the individual exposures taken on the 29'''. No 
fringe subtraction was done for the I band images. 

Surface photometry for KK98 250 was done using the 



Landolt (1983) \ 



ellipse fitting algorithm of Bender & Mollenhoff (1987) 



Ellipse fitting in the I band was restricted to the inner re- 
gions of the galaxy, since the outer low surface brightness 
emission is affected by fringing. The average ellipticity was 
found to be 0.72 ± 0.02 (corresponding to an inclination 
of ~ 79°, for an intrinsic thickness ratio go — 0.2). The 
position angle, measured between the north direction on 
the sky and the major axis of the receding half of the 
galaxy (see Sect. 13.3} . was found to be ^ 270°. The sur- 
face brightness profiles obtained from the ellipse fitting 
in V and I bands are well described by exponential disk 
with the scale length of 41.0" (-1.1 kpc) and 62.0" (-1.7 
kpc) respectively. The average <V-I> color is —1.2, with 
little variation with radius. The total magnitude after cor- 
recting for galactic extinction (using Ay^l.Ol mag and 
A/=0.59 mag; Schlegel et al. 1998) is 13.8±0.3 mag in V 
and 12.6±0.4 mag in I. No correction for internal extinc- 
tion was applied. 

KK98 251 has several bright stars superposed on it 
which could not be accurately removed. For this galaxy, 
we hence use the results of the I band photometry derived 
by Karachnetsev et al (2000)| who derived an I band ex- 
ponential scale length and a total I magnitude of 22.8" 
(—0.6 kpc) and 14.42 mag respectively. 

2.2. HI observations and analysis 

The GMRT (Swarup et al. 1991) observations of KK98 
250 and KK98 251 were conducted during the commis- 
sioning phase of the telescope. The field of view of the 
telescope (— 24 ) is large enough to cover both galaxies in 
a single pointing. The setup for the observations is given 
in Table n The calibrators were 3C48 and 3C286 (flux), 
2022-F616 (phase) and 3C286 (bandpass). 

The data were reduced in the usual way using standard 
tasks in classic AIPS. The GMRT has a hybrid configu- 



ration which simultaneously provides both high angular 
resolution (— 2 if one uses baselines between the arm an- 
tennas) as well as sensitivity to extended emission (from 
baselines between the antennas in the central array). Data 
cubes were therefore made using various (u,v) cutoffs, in- 
cluding 0-5 kA, 0-10 kA, 0-15 kA and 0-20 kA. The cor- 
responding angular resolutions are 43" x 38", 26" x 21", 
16" X 14" and 11.5" x 10". All the data cubes were decon- 
volved using the AIPS task IMAGR and then corrected 
for primary beam attenuation. 

Moment maps of the data cubes were made using the 
AIPS task MOMNT. Maps of the velocity field and line 
profile widths were also made in GIPSY, using single gaus- 
sian fits to the line profiles. The velocity fields produced 
by gaussian fitting were found to be in a reasonable agree- 
ment with that obtained from moment analysis. As we 
discuss in more detail later, we need to know the veloc- 
ity dispersion of the gas in order to estimate the pres- 
sure support in the HI disk. For dwarf galaxies, where 
the contribution from systematic rotation to the profile 
width in a single synthesized beam is small, the veloc- 
ity dispersion can generally be estimated from the width 
of the line profiles. However, for a highly inclined galaxy 
like KK98 250, the profile width is dominated by rota- 
tional motion even within a single synthesized beam. For 
this galaxy, we hence assume that the velocity dispersion 
is — 8 km s~^, a typical value for dwarf galaxies (e.g. 



Lake et al. (T990)l|Skillman et al. (1 988)). For KK98 251, 
from the widths obtained from gaussian fitting (which we 
prefer to use because the thresholding algorithm used by 
the MOMNT task results in an underestimation of the 
true profile width), we find that the profile widths vary 
from — 8.4 ± 0.7 km s^^ in the inner region (upto 80 ) to 
— 9.5 ± 0.9 km s~'^in the outer regions. 

Finally, a continuum image was made using the av- 
erage of the line free channels but no extended con- 
tinuum (to 3(7 limit of 1.8 mJy/beam for a beam size 
of 26 x 22 ) was detected from either galaxy. A high 
resolution (5.8 x 5.2 ) image also showed no compact 
sources associated with either galaxy down to a 3cr limit 
of 1.0 mJy/beam. 



Table 1. Parameters of the GMRT observations 



Parameters 



Value 



RA(2000) 20''30™33.0" 

Declination(2000) 60°21'17" 

Central velocity (heliocentric) 121.0 km s~^ 

Date of observations 22 June 2001 

Time on source 6 hrs 

Total bandwidth 1.0 MHz 

Number of channels 128 

Channel separation 1.65 km s~^ 

FWHM of synthesized beam 43" x 38", 26" x 21", 

16" X 14", 11.5" X 10" 
RMS noise per channel 3.5 mjy, 2.8 mjy, 2.2 mjy 

1.9 mJy 
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Fig. 1. [A] HI profile for KK98 250 obtained from 43 x 
38 data cube. The channel separation is 1.65 km s^^. 
Integration of profile gives a flux integral of 16.4 Jy 
km s-^and an HI mass of 12.1 x lO^M©. [B] HI pro- 
file for KK98 251 obtained from 43" x 38" data cube. The 
channel separation is 1.65 km s^^. Integration of profile 
gives a flux integral of 10.6 Jy km s^^and an HI mass of 
7.8 X IO^Mq. The dotted line shows a gaussian fit to the 
profile. 



3. Results and Discussion 

3.1. HI distribution 

The global HI emission profiles of KK98 250 and KK98 
251 (obtained from 43" x 38" data cube), are shown 
in Fig. H] The integrated flux of KK98 250 is 16.4 ± 
1.6 Jy km s~^, the (heliocentric) systemic velocity and 
the velocity width at 50% level of peak emission (W50) 
are 126.0±1.6 km s'^and 95.5 ± 2.3 km s'^ The sys- 
temic velocity is a good match to the single dish value 
of 127.0 ± 2.0 km s'^ HHuchtmeier et al. IMTJI . but the 
flux integral is considerably lower than the single dish 
value of 20.0 Jy km s-^ For KK98 251, the integrated 
flux is 10.6 ± 1.0 Jy km s^^, the systemic velocity is 
130.2 ± 1.7 km s^^and W50 is 51.7 ± 1.8 km s'^ Again, 
the systemic velocity is in reasonable agreement with the 
value of 126.0±2.0 km s~^obtained from single dish obser- 
vations (Hucht meier et al. 1997|l . but the integrated flux 
and W50 are considerably smaller than the single dish val- 
ues of 14.6 Jy km s~^and 64.0 km s~^ respectively. 

The GMRT fluxes could be lower than those obtained 
from single dish measurements either because of (i) a cal- 
ibration error or (ii) a large fraction of the HI being in 
an extended distribution that is resolved out. However, 
the flux of the point sources seen in the GMRT image 
are in good agreement with those listed in NVSS, indi- 
cating that our calibration is not at fault. Further, from 
our past experience in HI imaging of galaxies (with sizes 
similar to KK98 250 and KK98 251) with the GMRT, 
it seems unlikely that we have resolved out a large frac- 
tion of the total flux. Interestingly, a large discrepancy 
between the interferometric fluxes and the single dish 
fluxes was also seen in the DRAO images of these galaxies 
IjPisano fc Wilcots 2000|l . although the comparison in that 
case is complicated by the very poor signal to noise ratio of 



the DRAO data. Finally, we note that for both the galaxies 
there is a strong local HI emission at velocities very close 
to the systemic velocities. Hence, it is likely that the single 
dish integrated flux measurements were contaminated by 
blending of the HI emission from the galaxies with that 
of the galactic emission due to both the coarse velocity 
resolution ('^10.0 km s^^), as well as imperfect subtrac- 
tion of the foreground emission in the position switch- 
ing mode used in those observations. If we assume that 
the total fluxes are those measured at the GMRT, then 
the HI mass of KK98 250 is 12.1 ± 1.1 x 10^ Mq and 
Mni/Lg ~ 1.2. For KK98 251, the corresponding num- 
bers are 7.8 ± 0.7 x 10^ Mq and 1.6. 

Fig. El shows the integrated HI emission from KK98 
250 and KK98 251 at 26" x 21" resolution, overlayed on 
the digitized sky survey (DSS) image. Although KK98 250 
may be mildly warped, there is no clear signature of inter- 
action between the two galaxies. For KK98 251, two high 
density peaks are seen near the center. 

An estimate of the morphological center, position angle 
(PA) and inclination (assuming an intrinsic thickness ratio 
go — 0.2) of the galaxies were obtained by fitting elliptical 
annuh to the 43" x 38", 26" x 21" and 16" x 14" resolution 
integrated HI column density maps. For KK98 250, the 
inclination and PA were found to be 80 ± 4 and 267 ± 3 
degrees respectively. Due to the high density clumps in 
the center, ellipse fitting for KK98 251 was restricted to 
the outer contours of the 43" x 38" and 26" x 21" res- 
olution images. The inclination and PA of the HI disk 
were found to be 62 ± 5 and 220 ± 5 degrees respectively. 
For each galaxy, the value of the PA and inclination esti- 
mated from different resolution images match within the 
error bars. For KK98 250, the estimated parameters are 
also in good agreement with the values obtained from the 
optical image. 

Another input that we require for estimating the pres- 
sure support of the HI disk is the deprojected HI radial 
surface density profile. For highly inclined galaxies, de- 
projection using ellipse fitting does not lead to reliable 
estimates for the surface density. Hence, in the case of 
KK98 250, Lucy's (1974) iterative decovolution scheme, 
as adapted and developed by Warmels (1988)[ was used 
to derive the HI surface mass density. Fig.|2(AJ shows the 
best fit HI radial profile obtained from the 26" x 21" res- 
olution HI image using Lucy's scheme. Note that the de- 
convolved HI surface density profile shows a steep rise 
of ^40% in the central ~ 25 , i.e. within one synthe- 
sized beam. This is likely to be an artifact produced 
by the deconvolution; similar artifacts produced by this 
method have been seen earlier, by for e.g. Swaters (1999). 
We therefore fit the surface density Shi('') by a double 
Gaussian (shown as dotted fine in Fig.|3(A]), and use only 
the broad component for further analysis. For reference 
the full double Gaussian fit is given by: 



I]Hi(r) - Soi X e-^ /^^o^ + S, 



with Toi ~ 
5.0±0.1Mq pc 



73.9 ± 1.7 
"2 and E02 = 



02 X e 



ro2 



(1) 



14.4 ± 0.8 , E, 



01 



4.2 ± 0.2Mo pc- 
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Fig. 2. The optical DSS image of KK98 250 (north) and 
KK98 251 (south-east), shown in grayscales, with the 
GMRT 26 X 21 resolution integrated HI emission map 
(contours) overlayed. The contour levels are 0.02, 2.8, 5.7, 
8.7, 11.6, 14.4, 17.5, 20.1, 22.9, 25.8, 28.9, 31.7, 34.4, 37.4 
and 40.2 xl0^° atoms cm~^. 



For KK98 251, the HI surface density profile was ob- 
tained in the usual way from fitting elliptical annuli. The 
surface density Shi('') (shown in Fig. ElB]) was obtained 
from the integrated HI column density image at 26" x 21" 
resolution image and is well represented by a Gaussian: 



with ro = 34.2" ± 0.7" , c=19.2" ± 0.: 
O.IM0 pc-2. 



(2) 
and So == 7.8 ± 



Fig. 3. [A] The HI surface density profile of KK98 250 
derived from the HI distribution at 26 x 23 resolution 
(points). A double gaussian fit to the HI distribution is 
shown as dotted line. A decomposition of the double gaus- 
sian fit into a narrow and broad gaussian component is 
shown as dashed lines. See the text for more details. [B] 
The HI surface density profile of KK98 251 derived from 
the HI distribution at 26 x 23 resolution (points). The 
adopted gaussian fit is shown as a dotted line. 
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Fig. 4. The HI velocity fields of KK98 250 and KK98 
251 at 26 X 21 resolution. The contours are in steps 
of 5 km s~^and range from 80.0 km s~^to 170.0 km s~^. 

3.2. HI Kinematics 

The velocity fields of KK98 250 and KK98 251, derived 
from the 26" x 21" resolution HI data cube are shown in 
Fig. ^ The velocity fields are regular and the isovelocity 
contours are approximately parallel, which is a signature 
of rigid body rotation. 

The velocity field of KK98 251 is slightly lopsided, the 
isovelocity contours in the north-eastern half of the galaxy 
are more curved than the south-western half. Some small 
scale kinks are also seen in the isovelocity contours in the 
outer regions of the galaxy, indicative of a weak warping 
of the outer disk. These effects are more prominent in the 
higher spatial resolution velocity fields (not shown), which 
also show a mild twist in the kinematical major axis in the 
outer region of the galaxy, consistent with warping. 

3.3. HI Rotation Curve 

The rotation curves of the galaxies were derived separately 
for each of the 26 x21 ,16 xl4 and 11.5 xlO velocity 
fields, using the usual tilted ring model (Warner et al. 
1973). For the edge on galaxy KK98 250, as described in 
more detail below, the rotation curve was also derived by 
fitting to the position-velocity diagram. 

For each galaxy, the center and systemic velocity ob- 
tained from a global fit to the various resolution velocity 
fields matched within the error bars; the systemic velocity 
also matched with the values obtained from the global HI 
profiles. Keeping the center and systemic velocity fixed, we 
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fitted for the inclination and position angle (PA) in each 
ring. For KK98 250, keeping inclination as a free parame- 
ter in the tilted ring fit gave unphysical results, hence, the 
kinematical inclination of the galaxy was fixed to the value 
estimated from the HI morphology, viz 80°. Note that at 
such high inclination angles, the uncertainty in inclina- 
tion has only a small effect on the derived rotation curve 
(Begeman 1989). For example, if we fix the inclination to 
70°, the derived rotation curve is the same, within the er- 
ror bars. Fits for the position angle gave a value of r^ 267° 
(in good agreement with that derived from the optical im- 
age), with no systematic variation across the galaxy. The 
rotation curves (derived with the PA and inclination fixed 
to the values of 267° and 80°) are shown in Fig.E(A]. Note 
that the rotation curves derived from the different resolu- 
tion velocity fields match within the errorbars, suggesting 
that, in spite of being highly inclined, the effects of beam 
smearing are not significant. The final adopted rotation 
curve is shown in Fig. 13 A] as a solid line. 

For KK98 251, the inclination was found to be 65° 
(which agrees with that derived from ellipse fitting to the 
HI morphology, see Sect. I3.1|l . with no systematic vari- 
ation with the radius. The best fit position angle was 
^ 230° at all radii, except in the outermost regions of 
the galaxy where it changes to ^ 215°. Ignoring this rel- 
atively small change, (which in any case has negligible 
effect on the derived velocities), the rotation curve was 
derived by keeping the inclination and PA fixed at 65° 
and 230° respectively. Fig. ^A] shows the rotation curves 
for KK98 251, derived from the various resolution velocity 
fields. Note that all the rotation curves match within the 
errorbars. Recall that KK98 251 shows slight kinemati- 
cal lopsidedness. Consistent with this, the rotation curve 
derived independently for the approaching and receding 
sides are slightly different, although the difference is ev- 
erywhere less than 3.5 km s^^. For the purpose of mass 
modeling, we have used a mean of the rotation curves for 
the two sides. The errorbars on the mean rotation curve 
were obtained by adding quadratically the uncertainty re- 
ported by the tilted-ring fit as well as the difference in 
rotation velocities between the approaching and receding 
side. 

For highly inclined galaxies, rotation velocities derived 
from the titled ring fits to the velocity field could under- 
estimate the true rotation velocities, hence in such cases 
the rotation curve is often estimated by fitting to the high 
velocity edge of the emission (e.g. Sancisi & Allen 1979). 
While this method is well suited to large galaxies with 
flat rotation curves, it is not appropriate for galaxies with 
solid body rotation. Both the velocity field of KK98 250 
as well as the rotation curve derived from the tilted ring 
fit indicate solid body rotation, hence, instead of using the 
'edge fitting' technique, the rotation curve for KK98 250 
was derived by interactively fitting to the PV diagram, us- 
ing the task INSPECTOR in GIPSY (see FigElA]). The 
PA and the inclination in the interactive fitting were fixed 
to the values used in the tilted ring fit. The derived ro- 
tation curve matches, within the errorbars, to that de- 
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Fig. 5. [A] The rotation curves for KK98 250 derived 
from the intensity weighted velocity field at various reso- 
lutions. Circles, squares and triangles represent the rota- 
tion velocity derived from the 26 x 21 , 16 x 14 and 
11.5 X 10 resolution respectively. The adopted rotation 
curve is shown by solid line. [B] The adopted HI rotation 
curve for KK98 250 (dots) along with the Ha rotation ro- 
tation curve derived by de Blok et al. (2001) (triangles). 












.0..0 


B ..^ - 




30.00 


J' 




./f'-' 


».o 




10,00 


- /■■ : 




# 




■ t 


0,00 


- ,.,.," 



Fig. 6. [A] The rotation curves for KK98 251 derived 
from the intensity weighted velocity field at various resolu- 
tions. Circles, squares and triangles show the rotation ve- 
locity derived from the 26 x21 ,16 xl4 and 11.5 xlO 
resolution respectively. The adopted mean rotation curve 
is shown by dash-dot line. [B] The adopted mean rotation 
curve for KK98 251 (squares) and the rotation curve after 
applying "asymmetric drift" correction (circles). See the 
text for more details. 



rived from the tilted ring fit. As a further check of the 
robustness of the derived rotation curve, a model data 
cube for KK98 250 was constructed using the adopted ro- 
tation curve and the observed HI column density profile, 
with the task GALMOD in GIPSY. The model data cube 
was smoothed to a beam of 26" x 21" resolution and the 
moment maps were derived in the same manner as for 
the real data. Fig. ElB] shows the derived model velocity 
field for KK98 250. A residual (data- model) velocity field 
is shown in Fig. EIC] ; as can be seen, the model velocity 
field provides a good match to the observed field. 

Although KK98 251 is less inclined, a similar exercise 
of estimating the rotation velocities from the PV diagram 
was repeated for it. Fig|7[B] shows the adopted mean ro- 
tation curve projected onto the PV diagram - as can be 
seen, the mean rotation curve provides a reasonably good 
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Fig. 7. [A] The PV diagram of the galaxy along the kine- 
matical major axis for KK98 250, with the adopted rota- 
tion curve overlayed. [B] The PV diagram of the galaxy 
along the kinematical major axis for KK98 251, with the 
adopted rotation curve overlayed. 

fit to the data. A model data cube for KK98 251 was also 
constructed using the derived rotation curve, in the same 
manner as for KK98 250. Again, a good match between 
the model (not shown) and the observed field was found. 

The sensitivity of HI observations to the inner slope of 
the rotation curve has been the subject of much recent dis- 
cussion (e.g. van den Bosch & Swaters, 2001). One possible 
way of overcoming the relatively poor resolution offered 
by HI observations is to instead use Ha observations. For 
KK98 250, an Ha based rotation curve has been derived 
by Ide Blok et al. (2001) and is shown with triangles in 
Fig.[5[B]. As can be seen, although the Ha curve is steeper 
than the HI curve at intermediate radii, in the innermost 
regions of the galaxy, (where the effects of beam smearing 
are expected to be most severe), the two rotation curves 
in fact show an excellent agreement. 

As mentioned earlier, if one wishes to use the rotation 
curves to estimate the total dynamical mass, then one 
needs to account for the pressure support of the HI disk; 
this correction (generally called the "asymmetric drift" 
correction) is given by: 



V? = v^ - r X cr 



^[-(InSHi) + ^.{Ina') - :^(ln2h,)] , (3) 



^dr' 



dr 



dr' 



where Vc is the corrected circular velocity, Vq is the 
observed rotation velocity, a is the velocity dispersion, 
and hz is the scale height of the disk. Strictly speaking, 
"asymmetric drift" corrections are applicable to collision- 
less stellar systems for which the magnitude of the random 
motions is much smaller than that of the rotation veloc- 
ity. However, it is often used even for gaseous disks, where 
the assumption being made is that the pressure support 
can be approximated as the gas density times the square 
of the random velocity. The observed profile width can 
be used as an estimator of the velocity dispersion, after 
correcting for instrumental broadening and the finite size 
of the synthesized beam. This procedure doesn't work for 
very inchned galaxies, so, as discussed in Sect.|21for KK98 
250 we assume a constant velocity dispersion of 8 km s^^. 
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Fig. 8. [A]The observed velocity field of KK98 250 at 
26 X 21 resolution. [B]The model velocity field derived 
from the rotation curve at 26" x 21" resolution. The con- 
tours are in steps of 5 km s^^and range from 81.0 km s^^to 
171.0 km s~^[C] The residual velocity field ([A]-[B]). 
The contours are in steps of 1.0 km s~^and range from 
-4.0 km s~^to 4.0 km s"^. 



For KK98 251 we find that, after putting in these correc- 
tions, the estimated velocity dispersion is also ~ 8 km s~^. 
Further, in the absence of any measurement for hz, we as- 
sumed d(ln(hz))/dr = (i.e. that the scale height does 
not change with radius). Substituting these values back 
in the Eqn. |j2l and using the fitted Gaussian profile to 
the radial surface density distribution, the "asymmetric 
drift" correction was calculated and applied to the ob- 
served rotation velocities. For KK98 250 this correction is 
found to be small (less than 2.5 km s~^), compared to the 
errorbars on the rotation curve. On the other hand, the 
correction for KK98 251 is significant in the outer regions. 
The "asymmetric drift" corrected curve for KK98 251 is 
shown in Fig.ElB]. 
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Fig. 9. [A] Mass models for KK98 250 using the corrected rotation curve. The points are the observed data. The total 
mass of gaseous disk (dashed line) is 3.6 x lO^M0.The stellar disk (short dash dot hne) has Ty = 0.2, giving a stellar 
mass of 3.4 x lO^M©. The best fit total rotation curve for the constant density halo model is shown as a solid line, while 
the contribution of the halo itself is shown as a long dash dot line (the halo density is po — 37.6 ± 1.2 x 10"'^ M© pc~^ 
and core radius Yc = 1.5 ± 0.2). The best fit total rotation curve for an NFW type halo, using Ty = 0.0, c=1.0 and 
1^200=331.0 km s~^is shown as a dotted line. [B]The best fit mass model for an isothermal halo (solid line) to the 
hybrid rotation curve of KK98 250 consisting of the Ha rotation velocities in the center and the "asymmetric drift" 
corrected HI data in the outer regions of the galaxy. Also shown in the figure is the minimum disk fit mass model for 
an NFW type halo (dotted line). [C]Mass models for KK98 251 using the corrected rotation curve. The total mass of 
gaseous disk is 9.75 x IO^Mq. The stellar disk has T/ = 0.5, giving a stellar mass of 1.2 x IO^Mq. The best fit constant 
density halo model is po = 15.7 ± 1.3 x 10~^ M0 pc~^ and core radius Xc = 4.2 ± 1.6 (solid line). Also shown is the 



best fit total rotation curve for an NFW type halo, using T/ 
the text for more details. 



0.0, c=1.0 and U2oo=180.0 km s as a dotted line. See 



3 A. Mass Model 

In this section, we use the "asymmetric drift" corrected 
rotation curves, derived in the last section, to derive mass 
models for KK98 250 and KK98 251. 

The contribution of the stellar mass to the observed 
rotation curves were computed by assuming the galaxies 
to have exponential stellar disks, with a constant mass 
to light ratio (T) and an intrinsic thickness ratio (qo) 
of 0.2. We further assumed that the density distribu- 
tion in the vertical (2) direction falls off like sech^(2;/zo)i 
with zq independent of galacto-centric radius (see e.g. 
van der Kruit & Searle 1981, de Grijs & Peletier 1997). 
The contribution of the gaseous disks to the observed 
rotation curves were calculated using the observed HI 
surface mass density profiles, with the HI surface den- 
sity being scaled by a factor of 1.4 to account for the 
contribution from helium. There is little evidence that 
dwarf galaxies contain substantial amounts of molecular 
gas fe.g. Ilsrael et al. 19951 [Taylor et al. 1998| ), hence, no 
correction was made for molecular gas. We also neglected 
the contribution of ionized gas, if any. Since there is some 
evidence for similar vertical distributions of the HI and 
stellar disks (e.g. Bottema et al. 1986), we assumed the 
HI disks also to have a sech^(z/zo) vertical profile, with 
an intrinsic thickness ratio of go — 0-2. The circular ve- 
locities of the disk components were computed using the 



formulae given by Casertano (1983) For the dark matter 
halo, we considered two types of density profiles, viz. the 



modified isothermal profile and the NFW (Navarro et al. 
1996) profile. For mass models using a modified isother- 
mal halo, the free parameters are the halo central density 
Pq, core radius Tc and the mass to light ratio of the stellar 
disk, T. For the NFW models, the free parameters are the 
halo concentration parameter c, W200 (the circular velocity 
at the radius at which the halo density is 200 times the 
critical density) and the mass to light ratio of the stel- 
lar disk, T. Mass model were fit using the GIPSY task 
ROTMAS. 



3.5. KK98 250 

Since we could trace emission only from the brightest 
central regions of the galaxy in the I band, we use the 
more accurately determined V band scale length for the 
mass modeling. Fig. ElA] shows the best fit mass mod- 
els for KK98 250. The derived halo parameters are given 
in Table [3 For comparison, apart from the best fit mass 
model, the derived halo parameters are also given for the 
maximum disk, minimum disk and Ty — 0.7 (which was 
obtained from the observed color <V-I> of ~ 1.3 using 
the low metallicity Bruzual & Chariot SPS model using 
a modified Salpeter IMF; BeU & de Jong 2001). The to- 
tal dynamical mass of KK98 250, estimated from the last 
measured point of the rotation curve is 22.6 x lO^M© -at 
this radius more than 92% of the mass of KK98 250 is 
dark. 
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Table 2. Mass decomposition using isothermal halo 







KK98 250 






KK98 251 




HI Curve 


Tv 


re (kpc) 


po (10-3 M0 pc-3) 


Ti 


re (kpc) 


pa (10-3 Mo pc-3) 


Best fit 


0.20±0.05 


1.5 ±0.2 


37.6 ±1.2 


- 


- 


— 


Maximum disk 


3.0 


4.1 ±1.0 


10.0 ± 1.1 


1.0 


10.5 ±14.5 


13.4 ±1.6 


From observed <V-I> 


0.7 


1.6 ±0.1 


30.4 ± 1.0 


0.5 


4.2 ±1.6 


15.7 ±1.3 


Minimum disk 


0.0 


1.4 ±0.04 


39.0 ±1.0 


0.0 


2.8 ±0.5 


18.5 ±1.0 


Hybrid (HI+Hq) Curve 


0.2 


0.9 ±0.1 


65.4 ± 5.6 


- 


- 


— 


Hq Curve 


0.0 


0.63 ±0.08 


117.8 ±16.5 


- 


- 


- 



For mass models with an NFW halo, keeping Ty as 
a free parameter in the fit gave unphysical results. Even 
after setting Ty — 0, no reasonable fit could be obtained. 
As an illustration, Fig.|5lA] shows an NFW fit to the data, 
keeping the concentration parameters c fixed to 1, Ty=0.0 
and U200 chosen to minimize y^ . As can be seen, even at 
these extreme values for the parameters, the quality of fit 
is poor. 

We also fit mass models to a hybrid rotation 
curve (see Fig. EjB]), consisting of Ha data derived 
by de Blok et al. (200 1)[ in the inner regions of the galaxy 
and the "asymmetric drift" corrected HI rotation curve in 
the outer regions. Again, keeping Ty as a free parame- 
ter in the fit gave unphysical results, hence it was fixed 
to the value of 0.2, obtained from the best fit isothermal 
halo model, derived using the HI rotation curve alone. In 
any case, fixing Ty to a common value allows a mean- 
ingful comparison of the halo parameters derived using 
both the rotation curves. The derived halo parameters 
for the isothermal halo are given in Table El The ta- 
ble also shows the isothermal halo parameters derived 
by de Blok et al. (2001) using only the Ha rotation curve. 
We note that apart from the (probably not physically 
meaningful) maximum disk case, the halo parameters de- 
rived from the HI rotation curve are in good agreement 
with one another, but that they are substantially differ- 
ent from the parameters derived from the Hybrid or Ha 
rotation curves. In this context it is worth repeating, that 
the discrepancy between the Ha rotation curve and the HI 
rotation curve is largest at intermediate radii, and not at 
small radii as one would have expected, if the HI rotation 
curve suffered from beam smearing. 

Also shown in Fig.|n[B] is the Ty = 0, NFW halo fit to 
the Hybrid rotation curve. As can be seen, this fit overes- 
timates the observed rotation velocity in the inner regions 
of the galaxy; fixing Ty to more reasonable values worsens 
this discrepancy. The NFW halo parameters (correspond- 
ing to Ty = 0) are c=6.5±1.4 and U200 = 66.0 ± 13.0. 
The value of 6.5 for the concentration parameter is much 
lower than the range of NFW concentration parameters 
predicted by LCDM simulations (Bullock et al. 2001). It 
is worth mentioning here that even for the Ha rotation 
curve alone, isothermal halos provide a substantially bet- 
ter fit than an NFW halo (de Blok et al. 2001). 



3.6. KK98 251 

As discussed in section I^TI for KK98 251, we use the I 
band scale length derived by Karachentsev et al. (2000). 
For isothermal halo models, keeping T/ as a free param- 
eter in the fit gave negative values for T/. Further, if one 
keeps T/ fixed, the x^ continuously decreases as T/ is de- 
creased. Fig.EJC] shows the best fit mass model for KK98 
251 for a constant density halo using T/ = 0.5 (which 
corresponds to the observed <V-I> color, from the low 
metallicity Bruzual & Chariot SPS model using a mod- 
ified Salpeter IMF). The derived halo parameters using 
various values of T/ are given in Table |21 As can be seen 
from Table [21 the halo parameters are relatively insen- 
sitive to the assumed value of T/. The total dynamical 
mass (at the last measured point of the rotation curve) is 
Mt=9.5x1O*^M0. For KK98 251, we found that even for 
T/ = 0, with no value of c and i;2oo could a good fit using 
an NFW halo be obtained. The best fit NFW model is 
also shown in Fig. EtC] ; as can be seen, the data deviate 
substantially from the model. 



3.7. The Tully-Fisher relation 

It has been recently suggested that dwarf galaxies 
deviate systematically from the Tully-Fisher (TF) 
relation defined by bright galaxies (i.e. Stil (1999)1 
Swaters (1999)1 |McGaugh et al. (2000)| ), with small 
galaxies being underluminous compared to what would 
be expected had they followed the same TF rela- 
tion as ~ L* galaxies. Contrary to this suggestion. 



Pierini & Tuffs (1999) found no evidence of any break 
in the near-IR TF relation for faint dwarf galaxies. The 
reason for this discrepancy is unclear. We note however 
that high resolution HI images are probably crucial in 
studies of the TF relation for very faint galaxies. This 
is because (i) the inclination may often be difficult to 
obtain from images of faint irregular galaxies, and (ii) the 
50% HI velocity width may not be a good indicator 
of the rotation velocity in faint dwarf galaxies, where 
random motions are comparable to the peak rotational 
velocities (e.g. Camelopardalis B; [Begum et al. (2003)1 
DDO210; [ Begum fc Chengalur (2004)[ ). For such galax- 
ies, it is important to accurately correct for the pressure 
support ("asymmetric drift" correction), which is possible 
only if one has high resolution HI images. Unfortunately, 
there are only few dwarf galaxies with both I band pho- 
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tometry and HI synthesis imaging available. Interestingly 
enough, even without this correction, (or even a cor- 
rection for inclination), the HI velocity widths of dwarf 
galaxies with well measured distances do correlate with 
the absolute blue luminosity, albeit with large scatter 
(Huchtmeier et al. 2003). It is currently a matter of 
speculation as to whether the scatter would reduce 
significantly, once one puts in the corrections discussed 
above. 

We show in Fig. EH the TF relation for a sample of 
galaxies (with both HI synthesis imaging and I band pho- 
tometry) spanning a range of I magnitude from ^^ —23.0 
to —12.0. For the bright galaxies, the photometry and ro- 



tation velocities are from Verheijen (2001) while for the 
fainter galaxies, the I band magnitudes are from Patterson 
& Thuan (1996) and Makarova & Karachentsev (1998), 
while the rotation velocities were taken from Swaters 
(1999) and de Blok et al. (1996). The faintest two galax- 
ies (Camelopardalis B and DDO210) are from our ear- 
lier work (Begum et al. 2003; Begum & Chengalur 2004). 
For Camelopardalis B, the I magnitude is calculated us- 
ing the relation (B-V)t = 0.85(V-I)t-0.2 (Makarova & 
Karachentsev 1998). As can be seen from Fig. IIOL faint 
dwarf galaxies do tend to lie systematically below the 
TF relation defined by bright galaxies (dotted line). This 
trend persists even if one excludes the most discrepant 
bright galaxy (U7094, which has a poor quality rota- 



tion curve) in determining the bright galaxy TF relation 
(dashed curve). Note also that the scatter in the TF rela- 
tion increases at the faint end - this may be in part due 
to the uncertainty in the distances to the galaxies. 

4. Conclusions 

To conclude, we have presented high velocity resolution 
HI 21 cm synthesis images of the dwarf galaxies, KK98 
250 and KK98 251, as well as optical broad band images 
of KK98 250. We find that the HI disks of the galaxies do 
not show any clear signs of tidal disturbance and that both 
galaxies have regular velocity fields, consistent with rigid 
body rotation. We fit the rotation curves with dark matter 
halo mass models and find that both rotation curves can 
be fit using modified isothermal halos but not with NFW 
halos. Finally, for a small sample of galaxies with both 
I band photometry and HI synthesis images, we find, in 
agreement with earlier studies (which used single dish HI 
data) , that dwarf galaxies tend to lie below the I band TF 
relation defined by brighter galaxies. 
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